Disparity image fusion method for multiband stereo cameras

ABSTRACT

A disparity image fusion method for multiband stereo cameras belongs to the field of image processing and computer vision. The method obtains pixel disparity confidence information by using the intermediate output of binocular disparity estimation. The confidence information can be used to judge the disparity credibility of the position and assist disparity fusion. The confidence acquisition process makes full use of the intermediate output of calculation, and can be conveniently embedded into the traditional disparity estimation process, with high calculation efficiency and simple and easy operation. In the disparity image fusion method for multiband stereo cameras proposed by the method, the disparity diagrams participating in the fusion are obtained according to the binocular images of the corresponding bands, which makes full use of the information of each band and simultaneously avoiding introducing uncertainty and errors.

TECHNICAL FIELD

The present invention belongs to the field of image processing and computer vision, relates to multi-mode binocular disparity estimation and confidence calculation, and fusion and depth calculation of disparity information obtained by a multiband sensor, particularly to a confidence calculation method and a disparity image fusion method for multiband stereo cameras.

BACKGROUND

The binocular stereo vision technology based on visible light band is developed to be relatively mature. Visible light imaging has rich contrast, color and shape information, so the matching information between binocular images can be obtained accurately and quickly so as to obtain scenario depth information. However, visible light band imaging has defects, and the imaging quality thereof is greatly reduced, for example, in strong light, fog, rain, snow or night, which affects the matching precision. Therefore, the establishment of a multiband stereo vision system by using the complementarity of different band information sources is an effective way to realize three-dimensional reconstruction of space in special environments.

For example, a visible light band binocular camera and an infrared band binocular camera are used to constitute a multiband stereo vision system, and the advantage of not being affected by fog, rain, snow and light of infrared imaging is used to make up for the deficiency of visible light band imaging so as to obtain more complete and precise depth information.

In the multiband information complementary system, it is critical to fuse the multiband information to obtain the final depth information. At present, the research and application of a multiband stereo imaging system composed of infrared and visible light bands are more common. For example, Dr. Tong Ying fused infrared and visible light information at a two-dimensional image level, used infrared band information to repair and enhance the missing visible light band information, and conducted stereo matching and three-dimensional reconstruction on the basis of the generated fusion image (Research on Several Key Techniques of 3D Reconstruction Based on Infrared and Visible Light Image Fusion [D], Tianjin University, 2015). The key step of binocular depth estimation is to obtain the position, i.e., binocular matching information, of each pixel of the left figure L in the right figure R according to pixel intensity and relative position information between pixels, and the corresponding pixels of binocular images need to satisfy strict geometrical relationships so that the binocular depth information can be obtained according to the binocular matching information. Although the above method for fusing the visible light image and the infrared image at the two-dimensional image level recovers the approximate structure of the image and the pixel intensity information, the errors and uncertainty factors in the fusion process may change the geometrical relationships between matching pixel pairs. In addition, because the slight change of the pixel intensity and the slight change of the local structure caused by the uncertainty factors of the fusion may lead to the increase of the matching cost among the corresponding pixel points, the matching precision is influenced.

Another mainstream multiband information fusion mode is the fusion of point cloud information obtained in multiple bands, i.e., the infrared binocular image and the visible light binocular image as well as the calibration data of the infrared binocular camera and the calibration data of visible light are used for respectively generating point cloud corresponding to the infrared band and the visible light band and the fusion of the multiband binocular point cloud is explored by means of rigid transformation between an infrared unit and a visible light unit and various point cloud registration and fusion algorithms. (Research on Integration Algorithm of Infrared and Visible 3D Point Cloud). The above method respectively matches the binocular cameras in different bands, maps the images into three-dimensional information according to matching results, and fuse the information on a point cloud level, thereby greatly increasing the calculation amount and resource consumption and making it difficult to apply to actual products.

The present invention proposes a confidence calculation method and a multiband disparity information fusion method. The confidence information of each pixel is obtained in the stereo matching process of each band image; and after each band binocular camera respectively obtains disparity information with confidence, the information is directly fused on a two-dimensional disparity diagram according to the confidence and then the depth information is calculated according to the fused disparity diagram. The confidence information is obtained in the process of binocular disparity estimation. The disparity is directly obtained according to a single band binocular image, thereby avoiding the matching errors introduced by image fusion among the bands; and the confidence information is obtained from intermediate variables without adding extra calculation overhead. In addition, depth information fusion does not need the operation of conversion between two dimensions and three dimensions and projection mapping optimization, and greatly increases the efficiency compared with fusion on the point cloud information.

SUMMARY

The present invention aims to overcome the defects of the prior art and provides a disparity fusion method for the all-weather depth estimation process of microview and multiband cameras, which uses different band binocular cameras to obtain respective corresponding scenario disparities and fuses disparity diagrams obtained by sensors under different modes based on confidence strategy to obtain all-weather, more complete and more accurate depth information.

The technical solution of the present invention is as follows:

A disparity image fusion method for multiband stereo cameras comprises the following steps:

1. Calibrating respective internal and external parameters of multiband binocular cameras, wherein the internal parameters comprise a focal length and an optical center, and the external parameters comprise rotation and translation; correcting binocular images outputted by the binocular cameras in different bands into a parallel equipotential relationship.

2. Jointly calibrating the multiband binocular cameras to obtain position relationship parameters among devices of different bands, comprising rotation and translation.

3. Registering the multiband binocular cameras to obtain a coordinate conversion matrix of corresponding pixels among images collected by the devices of different bands, i.e., a homography matrix.

4. Matching the binocular images according to a semi-global matching (SGM) method to obtain respective cost diagrams of the multiband binocular images; and calculating the confidence and the disparity of each pixel point according to the cost diagrams to obtain disparity diagrams with the confidence.

5. Fusing the disparity diagrams obtained by cameras under two modes according to the confidence, rotation and translation parameters and the homography matrix to obtain more complete and more accurate depth diagrams; calculating a depth diagram corresponding to the fused disparity according to calibration parameters and the triangulation measurement principle.

The present invention has the following beneficial effects: the present invention proposes a strategy for fusing the disparity diagrams obtained by different band devices based on multiband binocular depth estimation, so as to achieve complementary advantages of different band devices. The present invention has the following advantages:

(1) The present invention proposes a method for obtaining pixel disparity confidence information by using the intermediate output of binocular disparity estimation. The confidence information can be used to judge the disparity credibility of the position and assist disparity fusion. The confidence acquisition process makes full use of the intermediate output of calculation, and can be conveniently embedded into the traditional disparity estimation process, with high calculation efficiency and simple and easy operation.

(2) The present invention proposes a disparity fusion method for multiband stereo cameras. The disparity diagrams participating in the fusion are obtained according to the binocular images of the corresponding bands, which makes full use of the information of each band and simultaneously avoiding introducing uncertainty and errors. The strategy of fusing the disparity according to the confidence can select the accurate information of each band and achieve the purpose of complementary advantages of multiband devices.

DESCRIPTION OF DRAWINGS

FIG. 1 is an overall flow chart of a solution;

FIG. 2 shows a detailed flow of a disparity estimation module;

FIG. 3 is a schematic diagram of influence of energy propagation on cost distribution;

FIG. 4 shows a detailed flow of a disparity fusion module; and

FIG. 5 is an effect diagram after disparity fusion.

DETAILED DESCRIPTION

The present invention fuses the disparity diagrams obtained based on a multiband device according to a multiband sensor device and a binocular disparity estimation method, and calculates distance information according to the triangulation measurement principle and the fused disparity to use the imaging advantages of devices of different bands under different environments. By taking the depth estimation of a pair of visible light binocular cameras and a pair of infrared binocular cameras as an example, a specific implementation solution is as follows:

FIG. 1 shows an overall flow of the solution.

I. The specific operation of a calibration and joint calibration module is as follows:

1. Respectively calibrating each lens of a visible light binocular camera and an infrared binocular camera and jointly calibrating the respective systems

1.1 Respectively calibrating the infrared camera and the visible light camera by the Zhangzhengyou calibration method to obtain internal parameters such as focal length and principal point position and external parameters such as rotation and translation of each camera.

1.2 Jointly calibrating the visible light binocular camera to obtain external parameters such as rotation and translation of two cameras of a binocular camera system; jointly calibrating the infrared binocular camera to obtain external parameters such as rotation and translation of two cameras of an infrared binocular system; correcting output image pairs according to the respective external parameters of the binocular camera system so that the binocular images outputted by the same binocular camera system satisfy the parallel equipotential relationship.

2. Jointly calibrating and registering the visible light binocular camera system and the infrared camera system

2.1 Jointly calibrating the left lens of the visible light binocular system and the left lens of the infrared binocular system by the Zhangzhengyou calibration method to obtain external parameters such as rotation and translation of the visible light camera and the infrared camera.

2.2 Simultaneously shooting the images of checkers in different planes by two pairs of binocular cameras, calculating the positional relationship of the same plane in the visible light image and the infrared image by using RT obtained by joint calibration and the detected checker corners, and representing the positional relationship with a homography matrix H.

II. The specific flow of a disparity estimation module is shown in FIG. 2. The specific operation is as follows:

1 Denoising and filtering input images by a Gaussian filtering algorithm. The weight of a Gaussian filtering window is determined by a Gaussian function (formula 1).

$\begin{matrix} {{h\left( {x,y} \right)} = e^{- \frac{x^{2} + y^{2}}{2\sigma^{2}}}} & (1) \end{matrix}$

wherein (x,y) is a point coordinate and a is a standard deviation. The Gaussian function is discretized to obtain a weight matrix, i.e., a Gaussian filter.

Through Gaussian filtering, noise can be effectively suppressed and the image is smoothed. Subsequent matching errors caused by the noise are prevented.

2. Matching the input images to obtain the disparity and the confidence. Stereo matching is mainly divided into the following processes: cost calculation, cost aggregation, confidence calculation and disparity calculation.

2.1 Firstly, conducting initial cost calculation. Sliding window matching based on Census features is taken as an example herein to describe the flow of matching cost calculation.

A Census feature descriptor of each pixel is obtained. A sliding window is used for search on a scanning line to calculate the cost corresponding to the possible disparity of each pixel (formula 2):

Cost_(d,d∈D) _(max) (x,y)=HD(CensusL(x,y),CensusR(x d,y))  (2)

In the formula, HD(⋅) represents a Hamming distance, and CensusL and CensusR are respectively Census feature descriptors of a left diagram pixel and a right diagram pixel. The output of cost calculation is a tensor of size height (H)×width (W)×maximum disparity (D), i.e., the cost diagram.

2.2 Conducting energy propagation on the initial cost diagram to obtain an energy diagram. The energy is described as follows:

$\begin{matrix} \left. {{E(D)} = {{\sum\limits_{p}{C\left( {p,D_{p}} \right)}} + {\sum\limits_{q \in N_{p}}{P_{1}{T\left\lbrack {{{D_{p} - D_{q}}} = 1} \right\rbrack}}} + {\sum\limits_{q \in N_{p}}{P_{2}{T\left\lbrack {{{D_{p} - D_{q}}} > 1} \right\rbrack}}}}} \right\rbrack & (3) \end{matrix}$

wherein C(p,D_(p)) is the cost at position p when the disparity is D_(p), T[⋅] is an indicator function, and the output is 1 when the input satisfies the conditions within [ ], otherwise, is 0. P₁ and P₂ are penalty terms. D_(q) is a disparity value at position q.

According to formula 4, in accordance with the global structure of the image, the cost distribution information of surrounding pixels is transmitted to the center pixel from multiple directions.

$\begin{matrix} {{L_{r}\left( {p,d} \right)} = {{C\left( {p,d} \right)} + {\min\ \begin{pmatrix} {{L_{r}\left( {{p - r},d} \right)},{{L_{r}\left( {{p - r},{d - 1}} \right)} + P_{1}},} \\ {{{L_{r}\left( {{p - r},{d + 1}} \right)} + P_{1}},} \\ {{\min\limits_{i}{L_{r}\left( {{p - r},i} \right)}} + P_{2}} \end{pmatrix}} - {\min\limits_{k}{L_{r}\left( {{p - r},k} \right)}}}} & (4) \end{matrix}$

L_(r)(p,d) is aggregation energy when the disparity at position p is d, and r is a transfer direction.

After energy propagation, a tensor of size height (H)×width (W)×maximum disparity (D), i.e., the energy diagram, is obtained. Energy propagation conducted successively in four directions of from top to bottom (TB), from left top to bottom right (LTB), from left to right (LR) and from right to left (RL) is taken as an example. The change of a cost curve of a pixel in the process of energy propagation is shown in FIG. 3.

2.3 Finding a disparity corresponding to minimum energy for each pixel position according to an energy propagation diagram, which is the integer disparity d(x,y) of the pixel.

$\begin{matrix} {{d\left( {x,y} \right)} = {\arg{\min\limits_{i \in D_{\max}}{{energy}(i)}}}} & (5) \end{matrix}$

The energy (⋅) is the energy after aggregation.

2.4 Calculating a subpixel level disparity diagram by using the energy diagram and an integral pixel disparity diagram.

III. Confidence calculation

In the calculation process of the binocular disparity, the confidence information of each pixel disparity can be calculated by the energy diagram information according to various confidence strategies and the confidence information is encoded into a binary string and connected in series with the disparity to save the storage space. The present invention proposes a confidence calculation strategy which comprises the following specific operation:

1 Searching the cost vector corresponding to each pixel in the energy diagram to obtain the minimum cost Cost_(min) corresponding to each pixel, corresponding disparity d_(min), second minimum cost Cost_(SecMin) and corresponding disparity d_(SecMmin);

2 Calculating the absolute value of a difference between two disparities:

abs2dis=|d _(min) −d _(SecMin)|  (6)

3 Calculating a ratio between two energy values:

Rate2Energy=Cost_(min)/Cost_(SecMin)  (7)

4 Judging a size relationship between two quantities. When abs2dis=thresh₁, if Rate2Energy<thresh₂, the confidence is set as 1, otherwise as 0; and when abs2dis>thresh₁, if Rate2Energy<thresh₁₃, the confidence is set as 1, otherwise as 0. thresh₁, thresh₂ and thresh₃ are obtained by statistics and stored in hardware in advance.

IV. Filtering the obtained disparity diagram to remove invalid disparity.

Firstly, speckle filtering is conducted on the image to remove the outlier in the disparity diagram.

$\begin{matrix} {{p\left( {x,y} \right)} = \left\{ \begin{matrix} {0,} & \begin{matrix} {\sum\limits_{{({i,j})} \in {\Omega{({x,y})}}}{T\left\lbrack {{p\left( {i,j} \right)} > {\left( {1 + t} \right) \cdot}} \right.}} \\ {{p\left( {x,y} \right)\left. {{p\left( {i,j} \right)} < {\left( {1 - t} \right) \cdot {p\left( {x,y} \right)}}} \right\rbrack} > {t\; 1}} \end{matrix} \\ {{p\left( {x,y} \right)},} & {other} \end{matrix} \right.} & (8) \end{matrix}$

p(x,y) is a disparity value at position (x,y); t and t1 are thresholds, counted and obtained by experiments and stored in hardware in advance; T[⋅] is an indicator function; and the output is 1 when the input satisfies the conditions within [ ], otherwise, is 0. Ω(x,y) is a local region centered on (x,y).

Median filtering is conducted on the disparity diagram. p(x,y)=median_((i,j)∈Ω(x,y))(p(i,j))  (9)

V. The specific operation of a disparity fusion module is as follows:

1. The disparity diagram disparity_(vis) obtained by the visible light binocular camera and the disparity diagram disparity_(ir) obtained by the infrared binocular camera are fused according to the homography matrix H, the translation and rotation positional relationships between the visible light system and the infrared system and two confidence marker bits. The specific operation is as follows:

1.1 Based on disparity_(vis), the coordinate of disparity_(ir) corresponding to each position is calculated according to the homography matrix, and a corresponding relationship between the coordinates of a pixel point N_(is) on disparity_(vis) and a pixel point p_(ir) matched with N_(is) on disparity_(ir) is as follows:

p _(vis) =H·p _(ir)  (10)

wherein H is the homography matrix formed by mapping the infrared disparity diagram to the visible light disparity diagram.

1.2 The disparity value and the confidence value of the matching point for the corresponding position are judged. If the disparity values of the corresponding positions of N_(is) and p_(ir) are both 0, no processing is conducted and disparity fusion is ended. If one of the disparity values is 0, the value of the corresponding position of disparity_(vis) is filled with non-zero party, and step 1.3 is conducted. If neither value is 0, the confidences of the corresponding positions of P_(vis) and p_(ir) are judged, and the corresponding position of disparity_(vis) is filled with the larger confidence. Step 1.3 is conducted.

1.3 Firstly, the disparity d_(ir) value corresponding to position p_(ir) (the coordinate is represented as (u_(ir), v_(ir)) under a corresponding pixel coordinate system) is transformed into a corresponding disparity value d_(vis) under a visible light binocular camera coordinate system, and then position P_(vis) (the coordinate is represented as (u_(vis), v_(vis)) under the corresponding pixel coordinate system) is filled. The coordinate at position P_(ir) under the corresponding camera coordinate system of the infrared device is represented as (X_(ir),Y_(ir),Z_(ir)).

p_(ir) is converted to the infrared camera coordinate system:

$\begin{matrix} \left\{ \begin{matrix} {X_{ir} = {\frac{Z_{ir}}{f_{ir}}\left( {u_{ir} - {cu_{ir}}} \right)}} \\ {Y_{ir} = {\frac{Z_{ir}}{f_{ir}}\left( {v_{ir} - {cv_{ir}}} \right)}} \\ {Z_{ir} = \frac{B_{ir} \cdot f_{ir}}{d_{ir}}} \end{matrix} \right. & (11) \end{matrix}$

wherein (cu_(ir),cv_(ir)) is the principal point position of a reference image of the infrared device. B_(ir) and f_(ir) are respectively the baseline length and the focal length of the infrared device. d_(ir) is a corresponding disparity value at position p_(ir).

According to the positional relationship between the visible light system and the infrared system, P_(ir) is converted to the visible light camera coordinate system:

[X′ _(vis) ,Y′ _(vis) ,Z′ _(vis)]^(T) =R·[X _(ir) ,Y _(ir) ,Z _(ir)]^(T) +t

(X′_(vis),Y′_(vis),Z′_(vis)) is the coordinate obtained by converting the point P_(ir) under the camera coordinate system of the infrared device to the camera coordinate system of the visible light device. R is the rotation relationship quantity between the infrared camera coordinate system and the visible light camera coordinate system, and t is the translation relationship quantity between the infrared camera coordinate system and the visible light camera coordinate system.

$\begin{matrix} \left\{ \begin{matrix} {X_{vis}^{\prime} = {\frac{Z_{u\;\_\;{vis}}^{\prime}}{f_{vis}}\left( {u_{vis} - {cu}_{vis}} \right)}} \\ {Y_{vis}^{\prime} = {\frac{Z_{v\;\_\;{vis}}^{\prime}}{f_{vis}}\left( {v_{vis} - {cv}_{vis}} \right)}} \\ {Z_{vis}^{\prime} = \frac{B_{vis} \cdot f_{vis}}{d_{z\;\_\;{vis}}}} \end{matrix} \right. & (12) \end{matrix}$

wherein (cu_(vis),cv_(vis)) is the principal point position of a reference image of the visible light device. B_(vis) and f_(vis) are respectively the baseline length and the focal length of the visible light device. d_(z_vis) is a disparity value obtained according to the coordinate Z′_(vis), Z′_(u_vis) is a depth obtained according to the coordinate X′_(vis), and Z′_(v_vis) is a depth obtained according to the coordinate Y′_(vis).

It can be obtained according to formula 12 that:

$\begin{matrix} \left\{ \begin{matrix} {Z_{u\;\_\;{vis}}^{\prime} = \frac{X_{vis}^{\prime} \cdot f_{vis}}{u_{vis} - {cu}_{vis}}} \\ {Z_{v\;\_\;{vis}}^{\prime} = \frac{Y_{vis}^{\prime} \cdot f_{vis}}{v_{vis} - {cv}_{vis}}} \\ {Z_{vis}^{\prime} = \frac{B_{vis} \cdot f_{vis}}{d_{z\;\_\;{vis}}}} \end{matrix} \right. & (13) \end{matrix}$

Thus, the disparity can be obtained according to formula 13:

$\begin{matrix} \left\{ \begin{matrix} {d_{u\;\_\;{vis}} = \frac{B_{vis} \cdot f_{vis} \cdot \left( {u_{vis} - {cu_{vis}}} \right)}{X_{vis}^{\prime} \cdot f_{vis}}} \\ {d_{v\;\_\;{vis}} = \frac{B_{vis} \cdot f_{vis} \cdot \left( {v_{vis} - {cv}_{vis}} \right)}{X_{vis}^{\prime} \cdot f_{vis}}} \\ {d_{v\;\_\;{vis}} = \frac{B_{vis} \cdot f_{vis}}{Z_{vis}^{\prime}}} \end{matrix} \right. & (14) \end{matrix}$

wherein d_(u_vis) is a disparity value obtained according to the X direction relationship, d_(v_vis) is a disparity value obtained according to the Y direction relationship, and d_(z_vis) is a disparity value obtained according to the Z direction relationship.

The disparity value that shall be filled at position v L vis is calculated.

d _(vis) =α·d _(u_vis) +β·d _(v_vis) +γ·d _(z_vis)  (15)

wherein α+β+γ=1 and α, β and γ are obtained through experiments.

d_(vis) is filled in position p_(vis).

Finally, the depth diagram is calculated according to the fused disparity, and the relationship formula between the disparity and the depth is as follows:

$\begin{matrix} {Z = \frac{B \times f}{d}} & (16) \end{matrix}$

wherein B is baseline length, f is the focal length, Z is the depth and d is the disparity. 

1. A disparity image fusion method for multiband stereo cameras, comprising the following steps:
 1. calibrating respective internal and external parameters of multiband binocular cameras, wherein the internal parameters comprise a focal length and an optical center, and the external parameters comprise rotation and translation; correcting binocular images outputted by the binocular cameras in different bands into a parallel equipotential relationship;
 2. jointly calibrating the multiband binocular cameras to obtain position relationship parameters among devices of different bands, comprising rotation and translation;
 3. registering the multiband binocular cameras to obtain a coordinate conversion matrix of corresponding pixels among images collected by the devices of different bands, i.e., a homography matrix;
 4. matching the binocular images according to a semi-global matching (SGM) method to obtain respective cost diagrams of the multiband binocular images; and calculating the confidence and the disparity of each pixel point according to the cost diagrams to obtain disparity diagrams with the confidence;
 5. fusing the disparity diagrams obtained by cameras under two modes according to the confidence, rotation and translation parameters and the homography matrix to obtain more complete and more accurate depth diagrams; calculating a depth diagram corresponding to the fused disparity according to calibration parameters and the triangulation measurement principle.
 2. The disparity image fusion method for multiband stereo cameras according to claim 2, wherein the specific process of the SGM method is as follows: 1.1 firstly, conducting initial cost calculation; when sliding window matching based on Census features is used, the process is: obtaining a Census feature descriptor of each pixel, and using a sliding window for search on a scanning line to calculate the cost corresponding to the possible disparity of each pixel: Cost_(d,d∈Dinax)(x,y)=HD(CensusL(x,y),CensusR(x d,y))  (2) in the formula, HD(⋅) represents a Hamming distance, and CensusL and CensusR are respectively Census feature descriptors of a left diagram pixel and a right diagram pixel; the output of cost calculation is a tensor of size height×width×maximum disparity, i.e., the cost diagram; 1.2 conducting energy propagation on an initial cost diagram to obtain an energy propagation diagram; the energy is described as follows: $\begin{matrix} \left. {{E(D)} = {{\sum\limits_{p}{C\left( {p,D_{p}} \right)}} + {\sum\limits_{q \in \mathcal{N}_{p}}{P_{1}{T\left\lbrack {{{D_{p} - D_{q}}} = 1} \right\rbrack}}} + {\sum\limits_{q \in \mathcal{N}_{p}}{P_{2}{T\left\lbrack {{{D_{p} - D_{q}}} > 1} \right\rbrack}}}}} \right\rbrack & (3) \end{matrix}$ wherein C(p, D_(p)) is the cost at position p when the disparity is D_(p), T[⋅] is an indicator function, and the output is 1 when the input satisfies the conditions within [ ]; otherwise, the output is 0; P₁ and P₂ are penalty terms; D_(q) is a disparity value at position q; according to formula (4), in accordance with the global structure of the image, the cost distribution information of surrounding pixels is transmitted to the center pixel from multiple directions; $\begin{matrix} {{L_{r}\left( {p,d} \right)} = {{C\left( {p,d} \right)} + {\min\begin{pmatrix} {{L_{r}\left( {{p - r},d} \right)},{{L_{r}\left( {{p - r},{d - 1}} \right)} + P_{1}},} \\ {{{L_{r}\left( {{p - r},{d + 1}} \right)} + P_{1}},} \\ {{\min\limits_{i}{L_{r}\left( {{p - r},i} \right)}} + P_{2}} \end{pmatrix}} - {\min\limits_{k}{L_{r}\left( {{p - r},k} \right)}}}} & (4) \end{matrix}$ L_(r)(p,d) is aggregation energy when the disparity at position p is d; r is a transfer direction; i and k respectively represent disparity values i and k at position p and position p−r; after energy propagation, a tensor of size height×width×maximum disparity, i.e., the energy propagation diagram, is obtained; 1.3 finding a disparity corresponding to minimum energy for each pixel position according to an energy propagation diagram, which is the integer disparity d(x,y) of the pixel; $\begin{matrix} {{d\left( {x,y} \right)} = {\arg{\min\limits_{i \in D_{\max}}{{energy}(i)}}}} & (5) \end{matrix}$ wherein the energy (⋅) is the energy after aggregation; 1.4 calculating a subpixel level disparity diagram by using the energy propagation diagram and an integral pixel disparity diagram.
 3. The disparity image fusion method for multiband stereo cameras according to claim 1, wherein the calculation method of the confidence specifically comprises: 2.1 searching the cost vector corresponding to each pixel in the energy propagation diagram to obtain the minimum cost Cost_(min) corresponding to each pixel, corresponding disparity d_(min), second minimum cost Cost_(SecMmin) and corresponding disparity d_(SecMmin); 2.2 calculating the absolute value of a difference between two disparities: abs2dis=|d _(min) −d _(SecMin)|  (6) 2.3 calculating a ratio between two energy values: Rate2Energy=Cost_(min)/Cost_(SecMin)  (7) 2.4 in judging a size relationship between abs2dis and Rate2Energy: when abs2dis=thresh_(ir) if Rate2Energy<thresh₂, the confidence is set as 1, otherwise as 0; when abs2dis>thresh_(ir) if Rate2Energy<thresh₁₃, the confidence is set as 1, otherwise as 0; wherein thresh₁ thresh₂ and thresh₃ are obtained by statistics and stored in hardware in advance.
 4. The disparity image fusion method for multiband stereo cameras according to claim 1, wherein the disparity fusion method is: the disparity diagram disparity_(A) obtained by a binocular system in band A and a disparity diagram disparity_(B) obtained by a binocular system in band B are fused according to the homography matrix H, the translation and rotation positional relationships between the binocular system in band A and the binocular system in band B and two confidence marker bits, which is specifically operated as follows: 3.1 based on disparity_(A), calculating the coordinate of disparity_(B) corresponding to each position according to the homography matrix, and a corresponding relationship between the coordinates of a pixel point p_(A) on disparity_(A) and a pixel point p_(B) matched with p_(A) on disparity_(B) is as follows: p _(A) =H·p _(B)  (10) wherein H is the homography matrix formed by mapping the disparity diagram obtained in band B to the disparity diagram obtained in band A; 3.2 judging the disparity value and the confidence value of the matching point for the corresponding position; if the disparity values of the corresponding positions of p_(A) and p_(B) are both 0, no processing is conducted and disparity fusion is ended; if one of the disparity values is 0, the value of the corresponding position of disparity_(A) is filled with non-zero party, and step 3.3 is conducted; if neither value is 0, the confidences of the corresponding positions of p_(A) and p_(B) are judged, and the corresponding position of disparity_(A) is filled with the larger confidence; step 3.3 is conducted; 3.3 firstly, transforming the disparity value d_(B) corresponding to position p_(B) into a corresponding disparity value d_(A) under a binocular camera coordinate system in band A, and then filling position p_(A); representing the coordinate at position P_(B) under the corresponding camera coordinate system of the system in band B as (X_(B),Y_(B),Z_(B)); wherein the coordinate under the corresponding pixel coordinate system at position p_(B) is represented as (u_(B),v_(B)), and the coordinate under the corresponding pixel coordinate system at position p_(A) is represented as (u_(A),v_(A)); converting p_(B) to the B camera coordinate system: $\begin{matrix} \left\{ \begin{matrix} {X_{B} = {\frac{Z_{B}}{fB}\left( {u_{B} - {cu_{B}}} \right)}} \\ {Y_{B} = {\frac{Z_{B}}{f_{B}}\left( {v_{B} - {cv_{B}}} \right)}} \\ {Z_{B} = \frac{b_{B} \cdot f_{B}}{d_{B}}} \end{matrix} \right. & (11) \end{matrix}$ wherein (cu_(B),cv_(B)) is the principal point position of a reference image of device B; b_(B) and f_(B) are respectively the baseline length and the focal length of the device B; Z_(B) is the coordinate in Z direction at P_(B); d_(B) is a disparity value corresponding to position p_(B); according to the positional relationship between the system in band A and the system in band B, converting P_(B) to the camera coordinate system of the system in band B: [X′ _(A) ,Y′ _(A) ,Z′ _(A)]^(T) =R·[X _(B) ,Y _(B) ,Z _(B)]^(T) +t (X′_(A),Y′_(A), Z′_(A)) is the coordinate obtained by converting the point P_(B) under the camera coordinate system of the device in band B to the camera coordinate system of the device in band A; R is a rotation relationship quantity between the camera coordinate systems of the device B and device A, and t is a rotation relationship quantity between the camera coordinate systems of the device B and the device A; $\begin{matrix} \left\{ \begin{matrix} {X_{A}^{\prime} = {\frac{Z_{u\;\_\; A}^{\prime}}{f_{A}}\left( {u_{A} - {cu_{A}}} \right)}} \\ {Y_{A}^{\prime} = {\frac{Z_{v\;\_\; A}^{\prime}}{f_{A}}\left( {v_{A} - {cv_{A}}} \right)}} \\ {Z_{A}^{\prime} = \frac{b_{A} \cdot f_{A}}{d_{z\;\_\; A}}} \end{matrix} \right. & (12) \end{matrix}$ wherein (cu_(A),cv_(A)) is the principal point position of a reference image of the device A; b_(A) and f_(A) are respectively the baseline length and the focal length of the device A; Z′_(u_A) is a depth value obtained from the coordinate X′_(A); Z′_(v_A) is a depth value obtained from the coordinate Y′_(A); d_(z_A) is a disparity value obtained according to the coordinate Z′_(A); from formula (12), obtaining: $\begin{matrix} \left\{ \begin{matrix} {Z_{u\;\_\; A}^{\prime} = \frac{X_{A}^{\prime} \cdot f_{A}}{u_{A} - {cu_{A}}}} \\ {Z_{u\;\_\; A}^{\prime} = \frac{Y_{A}^{\prime} \cdot f_{A}}{v_{A} - {cv}_{A}}} \\ {Z_{A}^{\prime} = \frac{b_{A} \cdot f_{A}}{d_{z\;\_\; A}}} \end{matrix} \right. & (13) \end{matrix}$ thus, from formula (13), obtaining the disparity: $\begin{matrix} \left\{ \begin{matrix} {d_{u\;\_\; A} = \frac{b_{A} \cdot f_{A} \cdot \left( {u_{A} - {cu_{A}}} \right)}{X_{A}^{\prime} \cdot f_{A}}} \\ {d_{v\;\_\; A} = \frac{b_{A} \cdot f_{A} \cdot \left( {v_{A} - {cv_{A}}} \right)}{Y_{A}^{\prime} \cdot f_{A}}} \\ {d_{z\;\_\; A} = \frac{b_{A} \cdot f_{A}}{Z_{A}^{\prime}}} \end{matrix} \right. & (14) \end{matrix}$ wherein d_(u_A) is a disparity value obtained according to the X direction relationship, d_(v_A) is a disparity value obtained according to the Y direction relationship, and d_(z_A) is a disparity value obtained according to the Z direction relationship; calculating the disparity value that shall be filled at position p_(A); d _(A) =α·d _(u_A) +β·d _(v_A) +γd _(z_A)  (15) wherein α+β+γ=1, and α, β and γ are obtained through experiments; filling d_(A) in position p_(A).
 5. The disparity image fusion method for multiband stereo cameras according to claim 3, wherein the disparity fusion method is: the disparity diagram disparity_(A) obtained by a binocular system in band A and a disparity diagram disparity_(B) obtained by a binocular system in band B are fused according to the homography matrix H, the translation and rotation positional relationships between the binocular system in band A and the binocular system in band B and two confidence marker bits, which is specifically operated as follows: 3.1 based on disparity_(A), calculating the coordinate of disparity_(B) corresponding to each position according to the homography matrix, and a corresponding relationship between the coordinates of a pixel point p_(A) on disparity_(A) and a pixel point p_(B) matched with p_(A) on disparity_(B) is as follows: p _(A) =H·p _(B)  (10) wherein H is the homography matrix formed by mapping the disparity diagram obtained in band B to the disparity diagram obtained in band A; 3.2 judging the disparity value and the confidence value of the matching point for the corresponding position; if the disparity values of the corresponding positions of p_(A) and p_(B) are both 0, no processing is conducted and disparity fusion is ended; if one of the disparity values is 0, the value of the corresponding position of disparity_(A) is filled with non-zero party, and step 3.3 is conducted; if neither value is 0, the confidences of the corresponding positions of p_(A) and p_(B) are judged, and the corresponding position of disparity_(A) is filled with the larger confidence; step 3.3 is conducted; 3.3 firstly, transforming the disparity value d_(B) corresponding to position p_(B) into a corresponding disparity value d_(A) under a binocular camera coordinate system in band A, and then filling position p_(A); representing the coordinate at position P_(B) under the corresponding camera coordinate system of the system in band B as (X_(B),Y_(B),Z_(B)), wherein the coordinate under the corresponding pixel coordinate system at position p_(B) is represented as (u_(B),v_(B)), and the coordinate under the corresponding pixel coordinate system at position p_(A) is represented as (u_(A),v_(A)); converting p_(B) to the B camera coordinate system: $\begin{matrix} \left\{ \begin{matrix} {X_{B} = {\frac{Z_{B}}{f_{B}}\left( {u_{B} - {cu_{B}}} \right)}} \\ {Y_{B} = {\frac{Z_{B}}{f_{B}}\left( {v_{B} - {cv_{B}}} \right)}} \\ {Z_{B} = \frac{b_{B} \cdot f_{B}}{d_{B}}} \end{matrix} \right. & (11) \end{matrix}$ wherein (cu_(B),cv_(B)) is the principal point position of a reference image of device B; b_(B) and f_(B) are respectively the baseline length and the focal length of the device B; Z_(B) is the coordinate in Z direction at P_(B); d_(B) is a disparity value corresponding to position p_(B); according to the positional relationship between the system in band A and the system in band B, converting P_(B) to the camera coordinate system of the system in band B: [X′ _(A) ,Y′ _(A) ,Z′ _(A)]^(T) =R·[X _(B) ,Y _(B) ,Z _(B)]^(T) +t (X′_(A), Y′_(A), Z′_(A)) is the coordinate obtained by converting the point P_(B) under the camera coordinate system of the device in band B to the camera coordinate system of the device in band A; R is a rotation relationship quantity between the camera coordinate systems of the device B and device A, and t is a rotation relationship quantity between the camera coordinate systems of the device B and the device A; $\begin{matrix} \left\{ \begin{matrix} {X_{A}^{\prime} = {\frac{Z_{u\;\_\; A}^{\prime}}{f_{A}}\left( {u_{A} - {cu_{A}}} \right)}} \\ {Y_{A}^{\prime} = {\frac{Z_{v\;\_\; A}^{\prime}}{f_{A}}\left( {v_{A} - {cv_{A}}} \right)}} \\ {Z_{A}^{\prime} = \frac{b_{A} \cdot f_{A}}{d_{z\;\_\; A}}} \end{matrix} \right. & (12) \end{matrix}$ wherein (cu_(A),cv_(A)) is the principal point position of a reference image of the device A; b_(A) and f_(A) are respectively the baseline length and the focal length of the device A; Z′_(u_A) is a depth value obtained from the coordinate X′_(A); Z′_(u_A) is a depth value obtained from the coordinate Y′_(A); d_(z_A) is a disparity value obtained according to the coordinate Z′_(A); from formula (12), obtaining: $\begin{matrix} \left\{ \begin{matrix} {Z_{u\;\_\; A}^{\prime} = \frac{X_{A}^{\prime} \cdot f_{A}}{u_{A} - {cu_{A}}}} \\ {Z_{v\;\_\; A}^{\prime} = \frac{Y_{A}^{\prime} \cdot f_{A}}{v_{A} - {cv}_{A}}} \\ {{Z_{A}\prime} = \frac{b_{A} \cdot f_{A}}{d_{z\;\_\; A}}} \end{matrix} \right. & (13) \end{matrix}$ thus, from formula (13), obtaining the disparity: $\begin{matrix} \left\{ \begin{matrix} {d_{u\;\_\; A} = \frac{b_{A} \cdot f_{A} \cdot \left( {u_{A} - {cu_{A}}} \right)}{X_{A}^{\prime} \cdot f_{A}}} \\ {d_{v\;\_\; A} = \frac{b_{A} \cdot f_{A} \cdot \left( {v_{A} - {cv_{A}}} \right)}{Y_{A}^{\prime} \cdot f_{A}}} \\ {d_{z\;\_\; A} = \frac{b_{A} \cdot f_{A}}{Z_{A}^{\prime}}} \end{matrix} \right. & (14) \end{matrix}$ wherein d_(u_A) is a disparity value obtained according to the X direction relationship, d_(v_A) is a disparity value obtained according to the Y direction relationship, and d_(z_A) is a disparity value obtained according to the Z direction relationship; calculating the disparity value that shall be filled at position p_(A); d _(A) =α·d _(u_A) +β·d _(v_A) +γd _(z_A)  (15) wherein α+β+γ=1 and α, β and γ are obtained through experiments; filling d_(A) in position p_(A). 